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We consider nonuniform black strings inside their event horizon. We present numerical evidence, 
that the singularity touches the horizon as the horizon topology changing transition is reached. 
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Introduction - Already Kaluza and Klein invoked a 
fifth dimension to unify gravity and electromagnetism, 
hiding the higher dimension subsequently by rendering 
it small and compact. Following these ideas, today su- 
pergravity and string theory also put forward the exis- 
tence of higher dimensions in their attempt to unify the 
fundamental forces of nature. 

In such _D-dimensional manifolds with p compact di- 
mensions, black holes can either be localized in the com- 
pact dimensions, or the black hole horizon can wrap the 
compact dimensions completely. Black holes which are 
localized in the compact dimensions have the horizon 
topology of a(D— 2)-sphere, S^^'^, and are called 'caged' 
black holes [HQ- contrast, when the horizon wraps 
the compact dimensions, the horizon topology reflects the 
topology of the compact manifold. In the simplest case, 
the single compact dimension is simply a circle, . The 
black holes then have the horizon topology of a torus, 
S^~^ X S^, and are referred to as black strings [U, 

As long as caged black holes are much smaller than 
the compact dimension, they resemble more or less D- 
dimensional Schwarzschild black holes (when they are 
static and uncharged). But as their size and with it 
their mass grows, they begin to feel the finite size of the 
compact dimension and deform accordingly, reaching a 
maximal size beyond which they no longer fit into the 
compact dimension [1, Black strings, on the other 
hand, exist for all values of the mass. However, as shown 
by Gregory and Laflamme (GL) such 'uniform' black 
strings (UBS), which do not depend on the compact co- 
ordinate, become unstable below a critical value of the 
mass. GL therefore suggested, that unstable UBS would 
decay to black holes, which possess higher entropy. 

In contrast, Horowitz and Maeda Q argued that the 
horizon of UBS could not pinch off in finite affine time. 
They therefore conjectured that the solutions would set- 
tle down to nontranslationally invariant solutions with 
the same horizon topology as the original configurations, 
'nonuniform' black strings (NUBS). Shortly after, such 
NUBS were found perturbatively [7] and numerically [1] . 
However, the NUBS cannot serve as endpoints of the in- 
stability, since they are too massive and possess too low 
an entropy 0, at least as long as the number of dimen- 
sions is smaller than 13 indicating that black holes 



are the end-state of the instability, nevertheless [l| . 

Emerging from the UBS branch at the marginally sta- 
ble solution, the NUBS form a branch of solutions along 
which the deformation of the horizon along the compact 
direction perpertually increases. This is illustrated in 
Fig. 1, where we exhibit the spatial embedding of the 
horizon of static D — 6 black string solutions for a se- 
quence of NUBS, yielding a geometrical view of the in- 
creasing nonuniformity of the solutions We here 
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parametrize the deformation by A — 2 \ 

where i?max and i?min represent the maximum and min- 
imum radius of a (Z? — 3)-sphere on the horizon. 
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FIG. 1: Spatial embedding of the horizon of D = 6 black 
strings with increasing deformation A, ranging from A = to 
A — 13.8. {Z denotes the proper length along the compact 
direction, R the proper radius of the horizon.) 



Extrapolating the numerical results to A — s- cx) shows 
that the maximal radius i?max assumes a finite limiting 
value, while the radius flmin of the 'waist' of the black 
strings shrinks to zero [4,|8|,|ia|. In the limit A ^ cx), the 
horizon will then pinch off, changing the horizon topology 



from 5^-3 X S'^ to 5^ 



^ at a singular 'merger' config- 
Study of the global and horizon 



uration [J, [O, d HI 
properties of NUBS and black holes in I? = 5 and D — 6 
dimensions indeed provides persuasive evidence that the 
NUBS branch merges with the black hole branch at such 
a topology changing transition [4, 10]. 

A deep question associated with the envisaged transi- 
tion is whether it is associated with the occurrence of a 
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naked singularity and thus a violation of cosmic censor- 
ship. Indeed, it appears possible, that as the black string 
pinches, a naked singularity is formed, because the sin- 
gularity, which orginally winds the compact dimension, 
gets broken pj]. To address this question, we here inves- 
tigate the inside of nonuniform black strings, to see how 
the curvature and in particular the singularity evolve, as 
A becomes large. Our results indicate, that the singu- 
larity will touch the horizon at the topology changing 
transition. 

Nonuniform Black Strings - We consider the Einstein 
action 
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for a D-dimensional spacetime with one compact direc- 
tion. We denote the periodic length of the compact di- 
rection by L. Here the last term is the Gibbons-Hawking 
surface term 

A convenient form of the metric is given by 



„2A . 



dr 
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-e^''f^dnl_, , (2) 



where z is the coordinate of the compact direction. A, 
B, and C are functions of f and z only, and / = 1 — 
{fo/f)^ '^- Thus the horizon resides at f = fp. The 
functions A, B, and C have to be determined from the 
Einstein equations Gl = 0, G^-hG^ = and G^ = 0. The 
remaining nontrivial equations G~ — G^ = and G!! — 
form the constraints. (For A = B — C ^ UBS arise.) 

We now focus on black strings in D = 6 dimensions, 
since they can be determined with great numerical accu- 
racy up to large deformations For solutions outside 
the horizon, f > fg, we introduce the coordinate f via 
f = ^/tq + f'^, while inside the horizon, f < fo, we define 
the coordinate r via f = v^^ — t^. In these coordinates 
the horizon resides at f = = r. Substitution in the line 
element Eq. ([2]) then yields for f < f o 

ds^ = e^\^df + e^^ i-dr^ + dz^) + e^^dill , (3) 

where e^^ = e'^^/{rl ~ r^), e^^ = e^^ , and e^^ = 
e^'^(fg — r^). Inside the horizon the Einstein equations 
are hyperbolic, the coordinate r playing the role of 'time'. 
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Initial conditions are given at the horizon (setting fp = 1) 
i(0, z) = An{z), B(0, z) = Bh(z), G(0, z) = Ch{z), 

drA\j.=o — dfA\r—Q — , drB\r=o — drB\fi—Q — , 
a,G|,=o = drC\r=0 = , (5) 

where the functions Ah(z), B}i{z) and Cuiz) are ob- 
tained from the solution outside the horizon. The origin 

of the z coordinate is chosen such that i?niin = e*^"'"' 
and i?,^,, = eC^H(L/2)^ 



Setting G* = G; G^ 



g: 



identities V^G'^ 



and V^G'^^ 
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G:^ = in the 



reveals that the 
constraints satisfy the advection equations 

{dr ± dz) {V^{Gl T (G; - Gl)/2)) = . 

Consequently, the constraints vanish everywhere, since 
they vanish at the horizon. 

Numerics - The hyperbolic differential equations (|4]) 
are re- written as a set of first order equations in r. For the 
evolution in r we use a fourth order Runge-Kutta method 
p^ . The 'spatial' coordinate is scaled via z z/L. The 
difference formulae for z derivatives are obtained from a 
polynomial approximation of the functions: If a function 
y has values yk at the gridpoints z^, we approximate the 
gth derivative of y at Zk by 



k-\-n 
j=k—n 
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Here pj^-* (z) is the gth derivative of the Lagrange poly- 
nomial of degree 2n, 
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To maintain the same order of approximation 
at all gridpoints we include auxiliary gridpoints 
{-2;„, -Zn-i, ■ • • , -Z2} and {I + zi,l + Z2, - ■ ■ ,1-1- z„}- 
The values of the functions at these points are ob- 
tained from the symmetry properties y{~z) — y{z) and 
y{l + z)^y{l-z). 

However, if a singularity is encountered at some point 
Zs, we have to restrict to the interval Zs < z < 1 — Zs- 
In this case the approximation of the z derivatives is less 
'symmetric'. Thus, if z^^ is the gridpoint next to z^, and 
ks < k < ks + n, we define 

Xq)(- \ N n Z - Zl 



3=ks 



l=k^ 
l¥^3 



Zj - Zl 



(7) 

and similarly for gridpoints close to 1 — Zs. 

As we argue below, the solution becomes singular if e*^ 
tends to zero. In order to avoid the singularity we have 
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to restrict to a domain, where e*^ is larger than some 
small number, typically « 10~^. Thus the computation 
is performed in two stages. In the first stage the solution 
evolves from the horizon (r = 0) until it becomes nearly 
singular at z = and z — I for some r = r^. During 
this stage we employ the difference formula ([6]). In the 
second stage the solution evolves from r = r^. In each 
Runge-Kutta step we check whether the solution is nearly 
singular at gridpoints Zk^ ■ If this is the case, we restrict 
to Zk^ < Zfe < 1 — Zk^ in the next Runge-Kutta step. 
During this stage the difference formula ^ is employed. 
The evolution stops when the z-interval shrinks to zero. 

Typical step sizes used are Ar « 1 — 4 x 10""* and 
the number of gridpoints ranges between N = 100 and 
N = 2400. The order of the difference formula is 2n = 
6 — 14. We have checked the consistency of the solutions 
for different choices of Ar, N and n, and also for different 
nonequidistant gridpoint distributions. 

Results - We have constructed NUBS solutions inside 
their horizon for several values of the deformation param- 
eter A. Near the horizon the solutions extend smoothly 
into the interior region. The function e'-^ decreases mono- 
tonically for fixed z. When tends to zero at some 
rs{z), the Kretschmann scalar diverges there, indicat- 
ing that the curvature singularity resides at rs{z). For 
UBS the radial coordinate of the singularity is constant, 
rs{z) = 1, while for NUBS with small deformation rs{z) 
shows an oscillation about r = 1 with small amplitude, 
which then increases with increasing A. We exhibit the 
coordinate z^ir) of the singularity in Fig. 2 for several 
values of the deformation parameter A. (Note the simi- 
larity with Fig. 1.) 
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FIG. 2: The coordinate Zs{r) of the location of the singularity 
is shown for several values of the deformation parameter A of 
black strings, ranging from A = to A = 13.8. 



As seen in the figure, the singularity approaches the 
horizon at z = and z = 1, when A becomes large. This 
is further demonstrated in Fig. 3, where we exhibit the 
minimal distance i?sing = J^^ e^^^'^^dr of the singularity 
from the horizon at the waist. We note that for large de- 
formation the distance decreases approximately inversely 
proportional to the deformation, i?sing '^-^ A~^. Thus we 



extrapolate that in the limit A oo the singularity will 
touch the horizon at the waist of the string, i.e., a naked 
singularity will appear at the topology changing transi- 
tion [i|,[i|. 




FIG. 3: The distance i?sing from the horizon to the singularity 
at the waist (z = 0) is shown versus the deformation A for a 
sequence of NUBS, together with its approximation at large 
A, -Rsing = A~^ (dotted). Also shown are the minimal horizon 
radius Rmin and its asymptote (dotted). 



We conclude that the deformation of the horizon re- 
flects itself in the location of the singularity. At the 
symmetric points z = and z — 1, where the function 
gCn yields the minimal horizon radius i?min, the distance 
i?sing of the singularity from the horizon is also minimal. 
In particular, for large A, Rmin ~ A~^ along with i?sing 
(see Fig. 3). Extrapolating to the limit of infinite defor- 
mation, we infer that the singularity touches the horizon 
at z = and z = 1, but is hidden well inside the horizon 
everywhere else. 

In order to get further insight into the geometry of 
space in the NUBS interior, we consider an isometric 
embedding of surfaces of constant Kretschmann scalar 
p.vpijR^^^'' ■ Thus, if such a surface has coordiantes 
rxiz), z, 9i, 02, f (at a fixed time), we define coordinates 
X„Zhy 



dXidX' + dZ^ 



„2B 



dR^ 



- {drK/dzf + 1 j dz^ + 
dZ'^ + R^dnl , (8) 



dnl 



where the Xi are expressed in terms of spherical coordi- 
nates in the second line. Regarding R and Z as functions 
of z then yields R{z) — e"-^ and 



Z{z) = L 



,2B 



1 — {drx/dz)' 



o2C 



dC/dz) dz' 



and 



where B and C are taken along the curve (r/f (z), z) 
the length scale L is reintroduced. 

We exhibit surfaces Z{R) of constant Kretschmann 
scalar K for NUBS with deformation A = 0.03 and 
A = 0.47 in Figs. 4 and 5, respectively [iB]. The limit 
of infinite K corresponds to the singularity at i? = 0. 
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FIG. 4: Isometric embedding of surfaces of constant 
Kretschmann scalar K for NUBS with A = 0.03. The straight 
lines correspond to surfaces of constant K of UBS with the 
same temperature. 
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FIG. 5: Same as Fig. 4 for A 
shown. 
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0.47. The horizon is also 



For UBS, which are also shown, R ~ K~'s, independent 
of Z. The NUBS surfaces, on the other hand, show an 
oscillation about constant R values, and thus a Z depen- 
dence of the power law. Here R tends more slowly to 
zero for Z = and faster for Z — ±Lk/'2, where Lk de- 
notes the periodic length of Z. Note, that Lk increases 
with increasing K . Clearly, the deformation of the hori- 
zon reflects itself also in the isometric embeddings of the 
surfaces of constant K. 

Outlook - Concerning the black hole - black string 
transition we have provided numerical evidence that at 
the merger point the NUBS singularity touches the hori- 
zon at the waist of the strings. What is still missing is the 
study of the interior of the caged black holes, as they ap- 
proach the merger point, to see how their curvature and, 
in particular, their singularity evolve. In D = 4 such a 
study has been performed recently , in higher dimen- 
sions, however, this remains a major numerical challenge. 

A GL instability arises also for rotating NUBS , and 
a similar instability occurs for asymptotically flat rotat- 
ing black holes in I? > 6 dimensions, which have a single 
angular momentum [l^. This instability then suggests, 
that as in the case of black strings, a branch of rotating 
'pinched' (nonuniform) black holes should arise at the 
marginally stable solution [l^. Exploring this analogy 
further leads to an intriguing phase diagram for rotating 
black holes, where horizon topology changing transitions 
from rotating pinched black holes to black rings, black 
saturns, and further configurations are expected to occur 
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